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Abstract 

Can supersymmetric models with a moderate stop mass be made consistent with the 
negative Higgs boson searches at LEP, while keeping perturbative unification manifest? The 
NMSSM achieves this rather easily, but only if extra matter multiplets filling complete SU (5) 
representations are present at intermediate energies. As a concrete example which makes use 
of this feature, we give an analytic description of the phenomenology of a constrained NMSSM 
close to a Peccci-Quinn symmetry point. The related pseudo-Goldstonc boson appears in 
decays of the Higgs bosons and possibly of the lightest neutralino, and itself decays into bb 
and TT. 



1 Introduction and motivations 



The absence of any clear signal of the Higgs boson(s) at LEP is a disturbing fact for the Minimal 
Supersymmetric Standard Model (MSSM). As well known a heavy stop could be the explanation. 
Although possible, however, this weakens the view that requires supersymmetry to be visible at 
the LHC, especially since the top, and the stop, have the strongest coupling to the Higgs boson 
system. In the attempt to avoid this quite unpleasant road, several proposals have in fact been 
made: among them, the consideration of the Next to Minimal Supersymmetric Standard Model 
(NMSSM) has received a great deal of attention. It is in fact true that the extra contribution to 
the quartic Higgs coupling arising in the NMSSM can easily accommodate a lightest Higgs boson 
even much heavier than in the MSSM[T],[2]. A strong constraint, however, to which we stick in 
this paper, is its compatibility with manifest perturbative unification. 

Although the NMSSM is a very minimal extension of the MSSM, it has a drawback: it allows 
to introduce several more parameters, which often make the various analyses difficult to follow or 
can even obscure the very search for significant phenomenological patterns. In this paper we try to 
clarify a possibility offered by the NMSSM to comply with the LEP constraint^ in a weakly fine- 
tuned and not too narrow region of its parameter space, while insisting on a relatively light stop. A 
key point is that such possibility rests on the largest possible values of the usual XSH1H2 coupling 
of the NMSSM consistent with manifest perturbative unification, including the possible existence 
of extra matter multiplets filling complete SU{5) representations at intermediate energies[8]. 

With CP conserved in the scalar sector, the NMSSM has three CP even and two CP odd 
neutral fields. With the standard definition of the Higgs doublets Hi and H2, the only scalar 
with tree level coupling to the vector boson pairs VV, often called h since it is the closest to the 
Standard Model (SM) Higgs boson, has the compositioiij^ 

h = h°cosp + h°smp (1.1) 

and tree-level diagonal mass-squared 

(m°)2 = M| cos^ 2(3 + \^v^ sin^ 2/3, (1.2) 

corrected by the well known radiative contribution from top-stop loops (for moderate mixing, 
< 1) 

-Ui<r^^.^o,'^.. (1.3) 



We do not show in this radiative correction a small positive contribution due to At, because in 
practice for moderate mixing it is compensated by a further negative 0{as<yt) correction at the 



two-loop order. Thus the one-loop result (1.3) remains a reasonable approximation. 



Note that equation (1.2) is valid for any scalar potential of the form 



V = V^™(i7i,ff2) + K{S)\Hi\' + ^ii{S)\H^\' - {fii{S)HiH2 + h.c.) + y\HiH2Y + V{S) (1.4) 



""^For other attempts see e.g. [3],|1],[5]. There are, however, two different aspects of the problem that may have 
not been equaUy addressed in these works: the level of fine-tuning in the Z-mass and the narrowness of the region 
of parameter space consistent with current data. 

^As usual, tan/3 = V2/V1, = {H°), vf + vj = 174 GeV. 
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Figure 1: Lower (red) curve: the maximal value of A(M^) as a function of tan (5 in NMSSM 
without extra matter at intermediate energies, subject to the condition AGUT/4vr < 0.3, 
fi^GUT = 0. Upper (blue) curves: same but with ns = 3 extra 5 + 5 at 1 TeV, and AGUT/4vr < 
0.3,0.15. 



i.e., in particular, for any NMSSM which stays perturbative up to the GUT scale. 

The phenomenology of the NMSSM in relation with the Higgs boson searches at LEP certainly 
depends on the value of m/^, but crucially also on the mixings of h with the two other CP even 
scalars, since being a diagonal entry of a positive definite squared mass matrix, gives only an 
upper bound on the mass squared of the lightest physical CP-even scalar. 

With this in mind, this paper consists of two logically independent but also complementary 
parts. In the first one we discuss the maximum possible values of the coupling A, and therefore of 
iTLh-, in presence of extra matter multiplets filling complete SU{5) representations at intermediate 
energies (Sect. [2]). Furthermore, based on the values that we find for rrih, we consider a simple and 
generic 2x2 mixing model between h and the lightest among the two remaining CP even scalars, 
which, before mixing with h, do not couple at all with VV (Sect. [s]). In the second part we 
describe a fully detailed and motivated version of the NMSSM with an approximate Peccei-Quinn 
sjTiimetry that realizes the phenomenological pattern outlined in the first part. This approximate 
symmetry restricts the number of effective parameters and makes possible an analytic description 
of most of the relevant features we want to underline. 



2 On the maximal value of the SH1H2 coupling 



From eq. (1.2) nih is especially sensitive to the value of the coupling A at the weak scale, which is 
constrained by demanding that A stays perturbative in its RGE evolution all the way up to the 
GUT scale. More specifically, since A grows with the energy from the weak to the GUT scale, we 
require for its value at the GUT scale, Xgut < 0.3 • 47r. 

The RGBs of A and of other relevant couplings can be found in Appendix |A} A significant 
indirect effect on the evolution of A is caused by the presence at intermediate energies of vector- 
like supermultiplets filling complete SU (5) representations [S]. These multiplets increase the gauge 
couplings at higher energies, which in turn slows down the growth of both A and yt, delaying the 
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Figure 2: The maximal value of rrih, see Eq. (1.3), in NMSSM without extra matter (red) 
and with = 3 extra 5 + 5 at 1 TeV (blue). The values of Aqut are the same as in Fig. [T| 
The stop mass is fixed at = 300 GeV with moderate mixing, \At/m^\ < 1. 



onset of nonperturbative behavior. This effect is illustrated in Fig. [T] which shows as function 
of tan f3 the maximum value of A at the weak scale without or with extra-matter effects (three 
(5 + 5) of SU{5) at the weak scale), for the current value of rrit = 171 GeVj^ Consequently, from 
eq.s ( 1.2[ 1.3[ ), Fig. [2] gives the maximum value of rrih for a moderate stop mass, = 300 GeV. 



The upper blue curves, for AGUT/47r = 0.3, 0.15 are again with three (5 + 5) of SU (5) at the Fermi 
scale, whereas the lower red curve, for AGUT/47r = 0.3, includes in the ROE evolution the standard 
matter effects only. 

Several features in these figures are worth being observed. All the curves in Fig. 2 go for large 
tan /3 to a common asymptotic value which is the upper bound on rrih in the MSSM (for the same 
rrif = 300 GeV and moderate mixing). Relative to this value, the increment in mh due to the 
extra three (5 + 5) is clearly significant, especially since without extra matter the maximum value 
of rrih barely touches the LEP bound on the SM Higgs boson mass of about 115 GeV. This is 
even more so since the upper limit on ruh is essentially saturated for wide variations of Aqut? in 
its upper range, as shown by the close upper curves in Fig. 2. 

Both Fig. 1 and 2 are for vanishing 15"^ coupling in the superpotential, but they are all 
insensitive to any choice of kqh^/Ait < 0.05, since n is rapidly driven to zero at lower energies by 
the RGE evolution. A larger kqut would however reduce the maximum A at the weak scale. 

Finally notice that the curves without any extra matter start at tan/3 ~ 1.6 because at lower 
tan f3, unlike in the case with extra matter, the top Yukava coupling hits by itself the perturbative 
bound of 0.3 ■ in before getting to the unification scale. 

Recently, Dine, Seiberg and Thomas [lOj have claimed that in a singlet extension of MSSM 
based on the superpotential ^HiH2 + \SHiH2 + \MS'^ one can raise the Higgs mass by a significant 
amount while maintaining manifest perturbative unification without extra matter at intermediate 
scales. E.g., one of their examples (^IQj, Sec. 4.1) had rrih = 120 GeV for tan/3 = 4, A = 0.7, 



■^The example 715 = 4 at 1 TeV emphasized in [8] gives non-perturbative values of acuT once the 2-loop terms 
are included in gauge beta functions. 
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715 = 


^5 = 3 




as{Mz) 


ac 


as{Mz) 


ac 




0.117 


0.041 


0.117 


0.103 


1-loop 


0.130 


0.043 


0.123 


0.154 


2-loop numerical 


0.129 


0.043 


0.122 


0.143 


2-loop analytical 



Table 1: Prediction for as{Mz) in the standard case {n^ = 0) and for 715 = 3 (5 + 5) at 1 
TeV. We use one and two loop gauge beta functions given in Appendix |A] without the two 
loop contributions of A and y, which can later be included perturbatively. The input MS 
values are d(M^)"i = 127.918, s'ln^ OwiMz) = 0.23122. We do not include any threshold 
corrections. The last line of the table is obtained by treating the 2-loop terms as perturbative 
corrections to the 1-loop results, following the standard method as described e.g. in |11] . 



rriz 



300 GeVQand no mixing, which is in clear contradiction with Fig. 2. As is stressed in the 



Introduction, our bound on rrih applies to any Higgs potential of the form (1.4), and in particular 
to the superpotential of ^0\. We believe that the expansion analysis of [lOj, based on integrating 
out S and analyzing the spectrum of light states in terms of coefficients of dimension 6 operators, 
must be breaking down, and this explains the discrepancy. 

We conclude this Section by analyzing the effect of extra SU{5) multiplets on the gauge 
coupling unification. In Table [l] we show the prediction of as{Mz) for ^5 = 3 from the running 
of the gauge couplings at one and two loops, compared with the standard case (ns = 0), without 
any threshold effect. In the same Table we give, for the two cases, the corresponding value of the 
unified coupling ac- As is well known, the one loop prediction is very close to the experimental 
value as{Mz) = 0.1176(20), and of course this conclusion is left unchanged by the addition of 
extra matter in full SU{5) multiplets. At two loops, the prediction for ^5 = 3 is brought closer 
to the experiment compared to the standard ris = result. However, the unavoidable presence of 
threshold corrections does not allow a significant distinction between the two cases. In fact, a i-th 
(5 + 5) split into a S'f/(3)-triplet of mass M^i and a SU{2) doublet of mass M^j, as{Mz) gives a 
further one loop threshold correction 

There is furthermore a two loop contribution from A itself, dominated by the UV, 

6as{Mz)^^ ^^s{Mz)^^JX^^^^Mc^\ ^2.2) 



as{Mz) 567r \ 27r^ 

i.e., numerically, 6as{Mz)/as{Mz) ~ -(1 ^ 2)% for AGUT/47r = 0.3 ^ 0.2. 
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Figure 3: For two reference values uih = 110 GeV (left) and nih = 120 GeV (right) we plot 



the normalized squared coupling, eq. (3.2 ), of the lightest scalar hi in the 2x2 mixing model 



eq. (3.1), plotted as a function of its mass mi and the mass of the heavier state m2- The 



region below the lower blue (upper red) curve is consistent with the 95% C.L. bounds [T6 
from nonobservation of hi at LEP2, assuming it decays into bb (bbbb). The heavier scalar /12 
to be consistent with the LEP2 searches, should have the mass above ~ 114 GeV. 



3 A simple 2x2 mixing model 

As already mentioned, what really matters for the NMSSM phenomenology, more than ruh itself, 
are the masses and compositions of the physical scalars. Before mixing with h, the MSSM has 
two other CP-even fields, si,S2 (in their mass-squared diagonal basis). Both their masses and 
compositions depend on all the various parameters of the NMSSM. Nevertheless, none of them is 
coupled to VV. 

Mixing of h with si and S2 (if one, or perhaps both, of these states are lighter than h) can help 
increase the mass of h. After the mixing, si and S2 acquire coupling to VV and become subject 
to LEP searches. As we are going to see, these mixing cannot be large for consistency with LEP. 
As such, one can analyze individually their additive effects without making any significant error. 
We can then consider a simplified 2x2 mixing modej^ between h and the lightest, si, among the 
two states not coupled to ZZ. Thus we consider a mass matrix 



-^2x2 



ml 

IS.VP? VP? 



2 



(3.1) 



with a fixed rrih and arbitrary m,,-^ < and Am^. 

In view of the previous Section and having in mind the LEP bound of about 115 GeV, only 
valid for the SM Higgs boson, we take for rrih two reference values, 110 and 120 GeV, close to 



^In fact they used the soft stop mass of 300 GeV, which corresponds to physical stop mass of 345 GeV, an 
irrelevant difference. 

^The effects of a 2x2 mixing model of this type has already been considered in the NMSSM [n],|S] and in the 
MSSM as well [12], [13], [H], [15]. 
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the upper bounds on rrih without or with extra matter respectively. In absence of mixing only 
the latter case would be compatible with LEP data. With mixing, however, which is generally 
present, the situation may change. 

In Fig.s 3, we describe the effect of mixing h with si in the two cases. In the plane of the 
two mass eigenvalues {mi, 1712) — from which we can uniquely reconstruct m^^ and Am^ — we give 
the isolines of the squared coupling of the lightest state to ZZ, normalized to the SM Higgs boson 
coupling: 

a.zz=f^y. (3.2) 
V 9hzz J 

From the data of Ref. [I6] this allows to determine in the same plane the 95% CL. bound from 
the non-observation of the lightest state, assumed to decay in bb with SM branching ratio. For 
later purposes we also consider the decay in bb bb with a branching ratio close to 1. Given the 
actual numbers, a quick way to understand from these figures the compatibility with LEP data is 
to see if there are values of the heaviest mass m2 above 115 GeV and simultaneously allowed by 
the bound on the lightest state. 

The conclusions are quite clear. With an unmixed value of rrih = 110 GeV, and a fortiori 
for lower values, it is hardly possible to obtain consistency with the LEP data |^ This means 
that, with a moderate stop mass and a small A^-term, the NMSSM without extra matter and with 
standard Higgs boson decays can perhaps be accommodated with LEP data, if at all, only in a 
small corner of its parameter space. This may explain the interest of considering the decay of the 
lightest state into rf rf, which is experimentally less constrained [3j. 

On the other hand, the ruh = 120 GeV case is obviously compatible with LEP data for small 
enough mixing. More important is that some mixing effects will inevitably be present, which can 
push the heavier state even further up with a somewhat reduced coupling to the ZZ, while keeping 
consistency with the LEP data for the lower state. This can be a characteristic feature of the 
NMSSM with extra matter contributing to the RGE running of the coupling constants, and is the 
phenomenological pattern to which we want to draw attention. 

4 An explicit example based on an approximate Peccei- 
Quinn Symmetry: PQ SUSY 

4.1 The Lagrangian and the allowed parameter space 

An independent motivation for the NMSSM is that it may provide a simple solution of the so called 
/x-problem: the supersymmetric superpotential mass term fiHiH2 gets replaced by X{S)HiH2 and 
all the mass terms in the Lagrangian originate from supersymmetry breaking. This possible 
solution of the /i-problem invites a symmetry explanation of the absence of mass terms in the 
superpotential. Such symmetries can be a continuous R-invariance and/or a Peccei-Quinn (PQ) 
symmetry. In this paper we choose a PQ symmetry since: i) it removes the ^S^ coupling, thereby 

^However, notice the point at (mi, 7712) ~ (95,115) GeV which has been emphasized in the hterature [13] in 
connection with a shght excess of events at LEP. 
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helping to maximize A at the weak scale; ii) it can reduce the number of parameters in the 
supersymmetry breaking Lagrangian as well, since PQ may be approximately realized in this 
sector without conflicting with experiments. This version of NMSSM, which we call "PQ SUSY," 
has a minimal number of parameters, and contains a light pseudo-Goldstone boson. For earlier 
considerations of NMSSM in the PQ limit, see |6],[7]. 

Up to the small breaking of the PQ symmetry, the Lagrangian is uniquely fixed by the super- 
potential term 

/ = \SH^H2, (4.1) 
by the soft non-supersymmetric piece of the scalar potential 



soft 



mg\S\ 



+ m\\Hi\^ + ml\H2\^ + {AxXSHiH2 + H.c.) 



(4.2) 



and by the gaugino mass terms, which we shall take large relative to X{S) (see below). Small 
breaking terms of the PQ symmetry, like 6V = m?S'^ + BfiHiH2+B..c. , will have to be present. 
However we assume them to be small enough only to give mass to the otherwise massless pseudo- 
Goldstone boson, without significantly affecting any of the remaining properties of the model. We 
have checked that this is a consistent approximation. 

When it exists, the CP-conserving, SU{2) x f/(l) U{1) breaking vacuum is related to the 
Lagrangian parameters by (x = m|/A^f^) 



Ml + 



l+x 



+ 



mr, 



cos 2(3 



(4.3) 



sin^ 2(3 = 2 



,m'( + m'2 + A^f^ 



l + x) - (l + x) 



(4.4) 



2A(1 



X 



sin 2/?. 



(4.5) 



Note that the scalar sector defined by eq.s ( 4.1|4.2 ) depends upon five parameters (apart from 



the pseudo-Goldtone mass mo)- A and m^, m|, Aa. We trade mj 2 for v and tan/9. A useful 
way to represent the various results, which we shall follow, is to show them in the plane (m^. Ax) 



for fixed values of A, tan/?. In particular the vacuum in eq.s (4.3, 4.4, 4.5) is indeed the true 
minimum of the overall potential only in a portion of this parameter space. Fig. |4] shows the 
allowed parameter space for tan/? = 1.5,2,2.5 and A close to the maximal allowed values from 
Fig. [TJ We see that Ax has a maximal and minimal allowed value for each ms in an interval 
< ms < m™*^^ ~ 70 GeV. Here m| > is required by the global stability. The upper limit on 
Ax comes from imposing that < at the minimum (4.3, 4.4, 4.5), so that it is preferred to the 
trivial stationary point at Vi = V2 = Vg = 0. The lower limit on Ax comes from the experimental 
bound on the chargino mass, m{x^) > 103 GeV [IT], via Eq. (4.14) below. The constraint of 
local stability does not further restrict the parameter space. 
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Figure 4: The allowed region of the {ms, Ax) plane (see the text) for tan (3 = 1.5, 2, 2.5 and 
A fixed at 0.65,0.7,0.75, respectively. 

4.2 Higgs boson and higgsino spectra 

The spectrum of the Higgs boson sector is straightforwardly obtained by expanding around the 
above minimum. For the single charged boson one finds 



1 + x 



\ 2 2 



(4.6) 



Out of the two neutral CP-odd states, one is massless in this approximation (the PQ Goldstone 
G) and the other has mass 

+AV(1 + x). 



1 + x 



Their compositions in terms of the neutral fields 



are 



G = — sin arts + cos a(cos (3^2 — sin /^tti), 
A = cos ails + sin a(cos Ptt2 — sin Pui), 



where 



tana 



A, 



(4.7) 
(4.8) 

(4.9) 
(4.10) 

(4.11) 



Xv{l + x) 

The 3x3 squared mass matrix of the CP-even neutral scalars is best written in the basis 

{H = cos /3/i2 — sin phi, h = cos phi + sin /3/i2, s) (4.12) 
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Figure 5: The masses of the two hghtest CP-even scalars in the model of Section |4j for 
(A, tan /3, j4a) = (0.7,2,400 GeV) and a range of nis values consistent with global stability 
of the scalar potential. We include the stop quantum correction with = 300 GeV. The 
region to the right of the dashed line {ms > 40 GeV) is excluded by LEP constraints on the 
S\ZZ coupling (see Fig. [6|. 



where it has the form 

+ (M| - \^v'^) sin^ 2/3 - |(M| - A2t;2) sin 4/3 -A^Xv cos 2/3 
A^^= [ * M| cos2 2/? + A V sin^ 2/3 -A^At; sin 2/3^^ | • (4.13) 

\ * * \^v^{\^x) 

Note, as anticipated, that ~ ("^h)^- Note also that one of the mixing terms between h 

and the two other scalars, A^i2 is always small, whereas the other, is essentially controlled 

by X (or m|). The mixing pattern discussed in the previous Section is precisely realized in this 
case, as shown in Fig. [5} where one has the two lightest scalar masses as functions of ms for 
(A, tan/?, y4;j) = (0.7,2,400 GeV). The lightest scalar mass is below the LEP limit; its dominant 
decay mode (see the next Section) is into 2 PQ pseudo-Goldstones: 5*1 — > GG 46. From Fig. 
|6] we see that for (A,tan/3) = (0.7,2) the LEP constraint on the S\ coupling to ZZ is satisfied 
in most of the parameter space allowed by the potential stability and the chargino mass bouncQ 
The heaviest CP-even scalar has mass ms-, 

^ Aa/(1 + x)i/2 ^i^jph for Ax = 400 GeV is in the 
380 -7- 400 GeV range, while A and are ~ 15 GeV heavier and lighter than 5*3, respectively. 

For heavy gaugino masses, the masses of the higgsinos are controlled by the effective fi- 
parameter 

/^ = Ay. = ^^^^sin2/?. (4.14) 
The single charged higgsino has mass m(x^) = /i, whereas the 3x3 neutralino mass matrix in 



^The processes e+e Z* ^ Si^2G followed by 6*1^2 GG could not possibly be seen at LEP2: the normalized 
squared couplings ^^5^ are tiny, < 10^^, one order of magnitude below the LEP2 limits [IB]. 
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A=0.7; tan/3=2 




ms [GeV] 

Figure 6: For (A,tan/3) = (0.7,2) we give tlie allowed regions in {ms,Ax) plane following 
from the stability of the potential and the chargino mass bomid (inside the blue curve, same 
as in Fig. |4]), and from the non-observation of Si GG 4b decays at LEP2 (to the left 
of the dashed red curve). In the same plot we show contours of the lightest CP-even scalar 
mass. As always, we include the stop quantum correction with = 300 GeV. 
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A=0.7, tan/3=2 
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Figure 7: The neutralino masses in tlie model of Section |4|for (A,tan/3) = (0.7,2) and the 
chargino mass fi in its typical range. We have Tn{xi) + "i(X2) = ^.(xs) since the neutralino 



mass matrix (4.15) is traceless. 



the basis {hi, h2, s) has the form 

(0 /i Afsin/3\ 
* Xvcosp] . (4.15) 
* * / 

Fig. [t] shows the values of the neutralino masses versus m(x^) in its typical range, 100 -i- 200GeV, 
for (A,tan/5) = (0.7,2). For these masses the LEP2 searches have not been possibly effective. 
Indeed, the process e"^e^ — > Z* ^ X1X2 is within the LEP2 kinematic limit in a part of the 
parameter space (see Fig. [T]). However, the production cross section turns out to be well below 
the ~ 0.1 pb limit set in [IT] due to phase space and coupling suppressions. 

Another possible process is e+e" — *• XiXi7isr with a photon (from Initial State Radiation) 
and missing transverse energy in the final state, which is constrained by LEP2 searches of extra 
neutrino species [IB]. However, we concluded that the existing data cannot rule out a Xi with a 
somewhat reduced ZxiXi coupling and mass above mz/2, as it is in our case (see also [IS]). 

If the gravitino is the lightest SUSY particle, the lightest neutralino will predominantly decay 
into the gravitino and the pseudo Goldstone boson G. However, if the SUSY breaking scale \/F 
exceeds about 1000 TeV, these decays happen outside the detector and do not modify collider 
phenomenology of the model (see Section [s]). 



4.3 Higgs boson couplings and branching ratios 

The phenomenology of the model is made peculiar by the presence of the light pseudoscalar G, 
with an unknown mass, rriG, coming from the breaking of the PQ symmetry and assumed to be 
relatively small. 

If niG is below the bb threshold, the pseudoscalar can be seen in radiative T decays [2D]- The 
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Production coupling 


Branching ratios 


^1 


^s^tt,^s,vv<20% (Fig. 8) 


BR{GG) > 98%) 


^2 


^s,u,^s,vv^ 100% 


See Figjgj 

BRixiXi) = 50 -^ 90% 
BR{GG) ^ 1 - BRixiXi) 


53 


^Sstt ^ 20%, ^53^^ negligible 


See Figjgj 

BRixiXj) - 35% (of which 50% into XiXi) 
BR{ZG) ~ 30% 
BR{SiSj) ^ 20% 



Table 2: The neutral CP-even Higgs boson dominant decay modes and the couphngs relevant 
for their production at the LHC via gluon fusion and vector boson fusion processes, for 
(A,tan/3,ylA) = (0.7,2,400 GeV) and < < 40 GeV. 



relevant branching ratio is given by [21 



BR{T 7G) = F(cosa)2(tan/3) 



2 Gprnl 



V^TTttEM 

2 X 10-^F(cosa)2(tan/3)^ F 



BR{T ^ 



mi 



where a is the angle in (4.11), and the suppression factor Fq < 0.5 is due to QCD, bound state 
and relativistic corrections (see [22], Section 3.1). For (A, tan /3, Aa) = (0.7,2,400 GeV) we get 
BR{T —>■ jG) ~ 0.5 X 10~^. The experimental limits on this branching ratio depends crucially on 
the decay properties of G. An interesting possibility occurs if 2mr < mc < 2mb, so that G decays 
into r+r^. In this case the current limit from CLEO is 



BR{T ^ 7G(^ tt)) < 10-^ 

and a dedicated run by BABAR may improve it soon by 1 — 2 orders of magnitude. 

Below we will assume that mc is above the bb threshold, corresponding to a relatively less 
restricted region of parameter space. The pseudoscalar then decays into bb and tt with branching 
ratios close to the branching ratios of the SM Higgs boson. 

All the couplings and decay rates for the other Higgs bosons are easily determined from the 
parameters of the model as given in the previous Section. Table |2] and Figs. [8]|9] illustrate the main 
features of the most relevant quantities for (A,tan/5, A^) = (0.7,2,400 GeV) and < mg < 40 
GeV (see Fig. g. 

Using these numbers one can make a preliminary conclusion that observing these states at the 
LHC will not be easy, since the production cross sections are suppressed, and the dominant decay 
products do not allow for easy background discrimination. Obviously, a more detailed study is 
required to assess the LHC discovery potential. 



12 



0.20 



0.15 



0.10 



0.05 



0.00 



A=0.7, tan/3-2, Aa=400 GeV 




Figure 8: The Sitt and SiVfj_V^ couplings squared of the lightest CP-even scalar 5i, nor- 



malized to the couplings of the SM Higgs boson analogously to Eq. (3.2). In this plot 
{X,tan(3,Ax) = (0.7, 2, 400GeV). 
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Figure 9: The dominant branching ratios of S2 (left) and 5*3 (right) for {X, tan f3, Ax) 
(0.7,2,400GeV). 
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4.4 Fine tuning 

It is interesting to know an estimate of the finetuning required to satisfy the various restrictions 
on the parameters of the model under consideration. For an early discussion of finetuning in the 
NMSSM see [23]. The first thing to check is if there is too strong dependence of the Z-mass, or 
of the vacuum expectation value v in eq. (4.3), on the various parameters. The strongest such 
dependence is on the parameter A\. Naively this would seem to require finetuning of the order 
Ax/X'^v'^ ~ 10 for Ax = 400 GeV. However, this estimate does not take into account the fact that 
the two cancelling terms are not totally independent: the variation of Ax influences the other term 



via the angle P as determined by the second equation (4.4). An estimate which takes this effect 
into account is given by the logarithmic derivativ^ 



A 



d log v"^ 



dhgA^ 



(4.16) 



TOj We see that tan/3 below ~ 1.7 starts to be 
= 2, Ax = 400 GeV is OK with less than 10% 



which can be evaluated numerically, see Fig. 
disfavored by this finetuning, although tan/5 
finetune. 

Whereas these considerations apply to the dependence of v on the low energy parameters, it is 
also necessary to check the consistency of the values of these same parameters with the expected 



contributions due to RGE evolutioEj^ Here the difference 
stop mass: 

o[m2 — m^) = ——[1717 — nif-) log — 



rrio 



gets a contribution due to the 



(no mixing j 



(4.17) 



TTii) I {ml — mf) = 1, see Fig 



11 



where log Amess/'i^ = 6 -h 40 for Amess = 100 TeV-i-10^^ GeV. This equation can be used to set an 
upper bound on the stop mass with no finetuning, i.e. for S{m: 
where we take low mediation scale, Amess = 
above, is safely within the allowed range. 
Finally we notice that, although v in eq 



100 Te\ 1" We see that m, = 300 GeV, as assumed 



(4.3) is only weakly dependent on rrig, the current 



experimental constraints mostly on the chargino mass (see eq. (4.14) and Fig. |4]) require ms to 
be below 40 50 GeV. At the same time there is a one-loop contribution to the running of m| 



due to Ay 



4tt 



^ \ 2 4 2 1 ^mess 

— A A;, log 



V 



(4.18) 



For A = 0.7, Ax = 400 GeV, and A = 100 TeV this gives 5m| = (110 GeV)^. Since the allowed 
range of m| is a factor 5 -h 10 smaller, to comply with this limit it is clear that the model under 
consideration would again prefer low mediation scale. 

®As a consequence of the above effect the derivative dv^/dA^ is actually negative. 

^Alternatively we could look directly at the dependence ofv on the high energy paremeters, which are considered 
more fundamental. 

^''To make this plot, parameters ml and m| have to be expressed in terms of A, tan /?, A\ and ms from equations 
(pl), (p|. 
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Figure 10: The logarithmic derivative Eq. (4.16) plotted as a function of Ax for (tan/3, A) 



(1.5, 0.65), (1.7, 0.7), (2, 0.7), (2.5, 0.75). Smaller tan (3 give bigger A (and hence require bigger 
finetuning). In this plot we assumed ms = 0, but the change for small allowed values of ms 
is negligible. 
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Figure 11: The maximal value of physical stop mass (with no mixing) which is consistent 
with the naturalness bound 6{m2 — m\)/{m'2 — mf) < 1, see Eq. (4.17). From below up: 



(tan/?. A) = (1.5, 0.65), (2, 0.7), (2. 5, 0.75). Mediation scale A^ess = 100 TeV. In this plot we 
assumed ms = 0, but the change for small allowed values of ms is negligible. 
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5 Conclusions and outlook 



Even if we assume, with good reasons indeed, that supersymmetry is relevant in nature, there is 
no water-tight argument that requires the presence of supcrsymmetric signals at the LHC. Our 
best hope is a natural solution of the hierarchy problem of the Fermi scale, which makes such a 
presence likely. For this to be the case, however, requires a low level of fine tuning in the Z-mass, 
i.e. a maximally natural solution of the hierarchy problem. The LEP limit on the Higgs boson 
mass is particularly important, since it excludes the most natural regions of parameter space of the 
simplest supcrsymmetric models. This amply motivates the focus on supcrsymmetric extensions 
of the SM that minimize this fine-tuning and remain, at the same time, reasonably simple. 

A particularly simple possibility for increasing the Higgs mass is to add a new quartic inter- 
action for the Higgs doublets via the superpotential interaction \SHiH2- However, to maintain 
perturbative unification of the gauge couplings, a clear success of weak scale supersymmetry, the 
value of A is limited. With minimal matter content this interaction provides at most an additional 
~ 10 GeV to the Higgs boson mass, leading to only a small allowed region of parameter space, 
even including mixing amongst the Higgs bosons. In contrast, with additional matter the pertur- 
bative evolution of couplings allows a larger value of A, increasing the Higgs boson mass by up to 
~ 20 GeV compared to the theory without the singlet field. Furthermore, in this case mixing can 
augment the Higgs boson mass by another 2-8 GeV, considerably enlarging the allowed region of 
parameter space. 

The extra matter implies that the gauge couplings arc larger in the UV, and the top coupling 
smaller, compared to the minimal matter case. Providing the extra matter fills complete SU(5) 
multiplets, the successful unification of gauge couplings at 1 loop is unaltered. The changes from 
2 loops and threshold corrections depends on the nature of the extra matter. In the case of 5 -|- 5 
representations, the prediction for as{Mz) from 2 loop running is decreased, improving the agree- 
ment with data, but this is offset by an expected increase in ag^Mz) from threshold corrections 
from non-degeneracies within the 5 + 5 multiplets. The situation with 10 + 10 representations is 
the opposite: 2 loop running increases the discrepancy with data, but is countered by the thresh- 
old corrections. In either case, the significance of gauge coupling unification is comparable to the 
standard case with minimal matter. 

In this work we have considered PQ SUSY — a version of the NMSSM that incorporates the 
above mechanism for enlarging the Higgs mass and is fully realistic, with a minimum number of 
parameters. The superpotential is assumed to be exactly invariant under a Peccei-Quinn symme- 
try, elegantly solving the problem, while the soft scalar interactions include small PQ breaking 
interactions to give a mass to the pseudo-Goldstone boson G. The theory possesses just two 
parameters more than the MSSM - one is the soft mass parameter for the scalar S, m|, and the 
other is the mass for G. A combination of vacuum stability and chargino mass limits implies a 
restricted range for m|, so that a 10 — 20% fine tuning is necessary, and a low messenger scale is 
preferred. 

The Higgs boson system has a few characteristic properties in its spectrum and in its couplings. 
The spectrum contains two CP-even neutral scalars relatively close in mass, one above and one 
below the "naive" LEP bound of 115 GeV by 10 ^ 20 GeV, and with a shared coupling to the 
vector boson pairs, VV. Only the sum of these coupling squared is close to the squared coupling of 
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the SM Higgs boson to VV. Related to the approximate Peccei-Quinn symmetry, the Higgs boson 
spectrum also contains a CP-odd light state, G, present in the main decay modes of all the CP 
even neutral scalars. G itself decays to hh and rf, with branching ratios close to those of a light 
SM Higgs boson. Quite clearly, to assess the discovery potential of such a Higgs boson system at 
the LHC or TeVatron requires, and deserves, a detailed examination. 

A small region of parameter space where 2m^ < mc < and G decays into r+r^ is also 
allowed; in this case the model predicts BR{T —>■ 7G) within reach of the existing B-factories. 

The general phenomenology of the model crucially depends on the properties of the lightest 
neutralino, which is predicted to have a mass near 50 GeV. Naturalness considerations suggest a 
low scale for supersymmetry breaking, \/F, so that the gravitino is the LSP and Xi the next-to- 
LSP. In this case, all superpartner production events at colliders will yield at least two Xiy with 
each decaying predominantly into a gravitino and a pseudo-Goldstone G, with a width of order 
10~^ eV for a/F = 100 TeV. Given the scaling T{xi) oc 1/-F^, Xi could therefore decay inside 
the detector if \/F is less than of order 1000 TeV. Pair production of xi LEP2 has a cross 
section of order 0.1 pb, leading in this case to events with 4 b jets and missing energy. We do 
not know if searches by the LEP experiments would have detected this signal. If not, the generic 
superpartner production signal at LHC/TeVatron may include 4 b jets together with the missing 
energy. The naturalness argument by itself is not sufficiently tight to prefer xi decays inside rather 
than outside the detector. 

Finally, the consistency of our model with the ElectroWeak Precision Tests merits further 
work; in particular, a non-negligible correction to the T parameter can be induced by values of A 
in the region of 0.6-0.8. 

It is important to note that there are alternative versions of the NMSSM with extra matter 
that incorporate both an enhanced Higgs boson mass, perturbative gauge coupling unification 
and a solution to the /i problem. The absence of mass parameters in the superpotential may be 
guaranteed by an R symmetry, that nevertheless allows the interaction {k,/3)S^ as well as XSH1H2. 
In order that the 5*^ interaction not substantially reduce the Higgs mass, the weak scale value of 
K should be less than about 0.1. However, the form of the renormalization group equations allows 
K at the unification scale to be close to unity, so this is not a powerful constraint on the theory. 
What is the form of the R symmetry breaking in the supersymmetry breaking scalar interactions? 
If Ax is the only significant R breaking parameter, and m| > 0, then this theory is a perturbation 
of the model discussed in this paper. On the other hand there is a new minimum for m| < that 
is very different from the one examined here, where the interaction prevents runaway behavior 
for Vs- There are also models of both the PQ and R types with large values of the symmetry 
breaking in the soft scalar interactions, but in these cases there are several more parameters that 
enter the phenomenology. Nevertheless, these models may be of interest since they may remove 
the need to tune m| to small values. 

Note added. After completion of this work we became aware of the work of P. Schuster and 
N. Toro [26] where the NMSSM in the PQ and in the R-symmetric limits is analyzed with special 
emphasis on the fine tuning issue. We believe that the present work usefully complements Ref. 
|26j in many different aspects. 
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A 2-loop beta functions 

For the convenience of the reader we give here 2-loop supersymmetric beta functions used to 
produce results in Section [2j The gauge coupling beta functions are: 



dgi 



dlogjj IGvr 



1 1 / 





f +^5 




■ 26 ■ 

5 




- 6 ■ 

5 




1+^5 


; bi-^top 


6 


, bi;\ = 


2 




-3 + ns 




4 








25 ^ 15 '^5 5 ^ S'^S 5 ^ 15 '^5 

I + 25 + Tns 24 
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5 + 9 14 + f n5 



Our result for b^j agrees with [25j. The dependence of bij on as given in jS] is wrong. 
The relevant beta functions of yt and of the NMSSM couplings A, k, are 



dyt ^ yt 
dlog/i IGtt^ 
d\ X 



dlogfi IGtt^ 

dK K 



{AX' + 3y^ + 2K'-gf-3gl) 



yt 



{22yf + 3X''yf + 3X^ + 2k' X^^ 



'IGtt 



2^2 



A 



(lOA^ + Qy^^A^ + Qyt + + 12X'k'^ 



dlogfi IGvr^ 



(6k^ + 6A^ 



K 



(167r 



2^2 



(167r2)2 

{24k^ + 2AX\' + 12A^ + ISy^X' ). 



In our analysis we have omitted the two loop contributions of gauge couplings to the running of 
ytiX, K. This is legitimate since gi do not approach non-perturbative values in the UV. 
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